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Abstract 


The strength of forward-backward (FB) multiplicity correlations is measured by the ALICE detector 
in proton-proton (pp) collisions at i/s = 0.9, 2.76 and 7 TeV. The measurement is performed in the 
central pseudorapidity region (| Tj | <0.8) for the transverse momentum GeV/c. Two separate 

pseudorapidity windows of width (Srj) ranging from 0.2 to 0.8 are chosen symmetrically around 
TJ = 0. The multiplicity correlation strength ((tcoit) is studied as a function of the pseudorapidity gap 
(Tjg^p) between the two windows as well as the width of these windows. The correlation strength is 
found to decrease with increasing Tjg^p and shows a non-linear increase with Srj. A sizable increase 
of the correlation strength with the collision energy, which cannot be explained exclusively by the 
increase of the mean multiplicity inside the windows, is observed. The correlation coefficient is also 
measured for multiplicities in different configurations of two azimuthal sectors selected within the 
symmetric FB Tj-windows. Two different conttibutions, the short-range (SR) and the long-range 
(LR), are observed. The energy dependence of bcorr is found to be weak for the SR component 
while it is sU'ong for the LR component. Moreover, the correlation coefficient is studied for particles 
belonging to various transverse momentum intervals chosen to have the same mean multiplicity. Both 
SR and LR contributions to bcorv are found to increase with pj in this case. Results are compared 
to PYTHIA and PHOJET event generators and to a string-based phenomenological model. The 
observed dependencies of bcorr add new constraints on phenomenological models. 
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1 Introduction 


We report a detailed study of correlations between multiplicities in pp collisions at 0.9, 2.76 and 7 TeV. 
The correlations are obtained from event-by-event multiplicity measurements in pseudorapidity (rj) and 
azimuth ((p) separated intervals. The intervals are selected one in the forward and another in the backward 
hemispheres in the center-of-mass system, therefore the correlations are referred to as forward-backward 
(FB) correlations. 

The FB correlation strength is characterized by the correlation coefficient, bcon, which is obtained from 
a linear regression analysis of the average multiplicity measured in the backward rapidity hemisphere 
((^b)«p) ^ function of the event multiplicity in the forward hemisphere (up): 

(«b)/1p =« + ^corr-«F ■ (1) 

This linear relation ([Til has been observed experimentally lll-fl] and is discussed in &0]. Under the 
assumption of linear correlation between Up and Ug, the Pearson correlation coefficient 


bcnrr — 


(«f) - 


( 2 ) 


can be used for the experimental determination of bcon 10 ] ■ Since the parameter a is given by a = 
(ug) — bcon{n^), it adds no additional information and usually is not considered 

Heretofore, FB multiplicity correlations were studied experimentally in a large number of collision sys¬ 


tems including e^e , ft^p, pp, pp and A-A interactions O.I4l8l- ll3h . No FB multiplicity correlations 
were observed in annihilation at y/s = 29 GeV. This was interpreted as the consequence of inde¬ 


pendent fragmentation of the forward and backward jets produced in this process lll4l . In contrast, in pp 
collisions at the ISR 113] at ^/s = 52.6 GeV 10] and in pp interactions at the SppS collider [15] sizeable 
positive FB multiplicity correlations have been observed. Their strength was found to increase strongly 
with collision energy yO, which was confirmed later at much higher energies {y/s > 1 TeV) in pp colli¬ 
sions by the E735 collaboration at the Tevatron il2ll and in pp collisions by the ATE AS experiment at 
the FHC {yfs = 0.9 and 7 TeV) lEl]. One of the observations reported by ATFAS is the decrease of ^corr 
with the increase of the minimum transverse momentum of charged particles. 

The STAR collaboration at RHIC analysed the FB multiplicity correlations in pp and Au-Au collisions 
at = 200 GeV ilTIl . Strong correlation was observed in case of Au-Au collisions, while in pp 
collisions bcon was found to be rather small (~0.1). In the present paper we relate this to the use of 
smaller pseudorapidity windows as compared to previous pp and pp measurements. 


Forward-backward multiplicity correlations in high energy pp and A-A collisions also raise a consider¬ 
able theoretical interest. First attempts to explain this phenomenon |l7,18-^ were made in the frame¬ 
work of the Dual Parton Model (DPM) | 0 ] and the Quark Gluon String Model (QGSM) |21, 22]. They 
provide a quantitative description of multiparticle production in soft processes. In improved versions 
of the models, collectivity effects arising due to the interactions between strings, which are particularly 
important in the case of A-A interactions, were taken into account ll23l426ll. These effects are based on 
the String Fusion Model (SFM) proposed in ll27l. l28ll. It was shown that these string interactions lead to 
a considerable modification of the FB correlation strength, along with the reduction of multiplicities, the 
increase of mean particle pt> and the enhancement of heavy flavour production in central A-A collisions 

ESQ. 


FB correlations are usually divided into short and long-range components 0 , 0 ]. In phenomenological 
models, short-range correlations (SRC) are assumed to be localized over a small range of T]-differences, 
up to one unit. They are induced by various short-range effects from single source fragmentation, in¬ 
cluding particles produced from decays of clusters or resonances, jet and mini-jet induced correlations. 
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Long-range correlations (LRC) extend over a wider range in rj. They originate from fluctuations in 
the number and properties of particle emitting sources (clusters, cut pomerons, strings, mini-jets etc.) 

0,0, Ellis. 


The SFM predicts that the variance of the number of particle-emitting sources (st rings) s hould be damped 
by their fusion, implying a reduction of multiplicity long-range correlations |23, 25, H]. Contrary to this 
prediction, long-range correlations arising in the Color Glass Condensate model (CGC) OlH have been 
shown to increase with the centrality of the collision 11321] . Therefore, the investigation of correlations 
between various observables, measured in two different, sufficiently separated T]-intervals, is considered 
to be a powerful tool for the exploration of the initial conditions of hadronic interactions 1^ . In the case 
of A-A collisions, these correlations induced across a wide range in rj are expected to reflect the earliest 
stages of the collisions, almost free from final state effects 0211341] . The reference for the analysis of 


A-A collision dynamics can be obtained in pp collisions by studying the dependence of FB correlations 
on collision energy, particle pseudorapidity, azimuth and transverse momenta. 


This paper is organized as follows: Section [2] provides experimental details, including the description of 
the procedures used for the event and track selection, the efficiency corrections and systematic uncer¬ 
tainties estimates. Sections [3] and |4] discuss the results on FB multiplicity correlation measurements in 
T] in pp collisions at y/s = 0.9, 2.76 and 7 TeV and in T]-0 windows at y/s = 0.9 and 7 TeV. In Section 
[3l we present dependences of the correlation coefficient on the gap between windows, their widths and 
the collision energy. In Section |4l multiplicity correlations in windows separated in pseudorapidity and 
azimuth are studied, and the comparison with Monte Carlo generators PYTHIA6 and PHOJET is dis¬ 
cussed. Results on multiplicity correlations in different pj ranges in pp collisions at y/s = 1 TeV are 
presented in Section 


2 Data Analysis 

2.1 Experimental setup, event and track selection 


The data presented in this paper were recorded with the ALICE detector ll35h in pp collisions at yfs = 0.9, 
2.76 and 7 TeV. Charged primary particles are reconstructed with the central barrel detectors combin¬ 
ing information from the Inner Tracking System (ITS) and the Time Projection Chamber (TPC). Both 
detectors are located inside the 0.5 T solenoidal field. 


The ITS is composed of 3 different types of coordinate-sensitive Si-detectors. It consists of 2 silicon 
pixel innermost layers (SPD), 2 silicon drift (SDD) and 2 silicon strip (SSD) outer detector layers. The 
design allows for two-particle separation in events with multiplicity up to 100 charged particles per cm^. 
The SPD detector covers the pseudorapidity ranges |t]| < 2 for inner and |t]| < 1.4 for outer layers, 
acceptances of SDD and SSD are |t] | < 0.9 and |t] | < 1, respectively. All ITS elements have a radiation 
length of about 1.1% Xq per layer. The ITS provides reliable charged particle tracking down to transverse 
momenta of 0.1 GeV/c, ideal for the study of low-px (soft) phenomena. 

The ALICE TPC is the main tracking detector of the central rapidity region. The TPC, together with the 
ITS, provides charged particle momentum measurement, particle identification and vertex determination 
with good momentum and dE/dx resolution as well as two-track separation of identified hadrons and 
leptons in the pj region below 10 GeV/c. The TPC has an acceptance of |? 7 1 < 0.9 for tracks which reach 
the outer radius of the TPC and up to |t]| < 1.5 for tracks that exit through the endcap of the TPC. 

For the present analysis, minimum bias pp events are used. The minimum-bias trigger required a hit in 
one of the forward scintillator counters (VZERO) or in one of the two SPD layers. The VZERO timing 
signal was used to reject beam-gas and beam-halo collisions. The primary vertex was reconstructed 
using the combined track information from the TPC and ITS, and only events with primary vertices 
lying within ±10 cm from the centre of the apparatus are selected. In this way a uniform acceptance in 
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Fig. 1: Illustration of the valuables drj, and Fig- 2: Illustration of sets of p-windows with different 
used in the present analysis. widths 5p and separation gaps Pgj,p. 
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Fig. 3: Illustration of 8 configurations of azimuthal sectors. Forward and backward pseudorapidity windows of the 
width 5p = 0.2 are additionally split into 8 azimuthal sectors with the width 5(p=7t/4. The red sectors correspond 
to the first window of the FB pair, the green sectors to the second one. The variable ^j,gp is the separation in 
azimuthal angle between centres of the sectors. 

the central pseudorapidity region |p| < 0.8 is ensured. The data samples for ^/s = 0.9, 2.76 and 7 TeV 
comprise 2 x 10^, 10 x 10^, and 6.5 x 10^ events, respectively. Only runs with low probability to produce 
several separate events per one bunch crossing (so-called pile-up events) were used in this analysis. 

To obtain high tracking efficiency and to reduce efficiency losses due to detector boundaries, tracks 
are selected with pj > 0.3GeV/c in the pseudorapidity range |p| < 0.8. Employing a Kalman filter 
technique, tracks are reconstructed using space-time points measured by the TPC. Tracks with at least 70 
space-points associated and track fitting X^/^dof less than 2 are accepted. Additionally, at least two hits 
in the ITS must be associated with the track. Tracks are also rejected if their distance of closest approach 
(DCA) to the reconstructed event vertex is larger than 0.3 cm in either the transverse or the longitudinal 
plane. Eor the chosen selection criteria, the tracking efficiency for charged particles with pj > 0.3 GeV/c 
is about 80%. 



2.2 Definition of counting windows 

Two intervals separated symmetrically around rj = 0 with variable width Sr] ranging from 0.2 to 0.8 
are defined as “forward” (E, rj > 0) and “backward” (B, T] < 0) . Correlations between multiplicities of 
charged particles (n) are studied as a function of the gap between the windows (denoted as figap). Another 
convenient variable is which is the separation in pseudorapidity between centres of the windows. 
These variables are illustrated in Eig. [T] and all configurations of window pairs chosen for the analysis 
are drawn in Eig. |2j 

The analysis is extended to correlations between separated regions in the rj-cp plane (sectors). The (p- 
angle space is split into 8 sectors with the width 5(p=7i/4 as shown in Eig. [3] This selection is motivated 
by a compromise between granularity and statistical uncertainty. The definitions and equations, described 
in Section [TJ remain the same for the rj-cp windows. The acceptance of the windows is determined by 
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Fig. 4: Illustrative example of forward versus backward raw multiplicity distribution for windows with 5rj — 0.6 
and Tjg^p = 0.4 at ^ = 1 TeV (left) and corresponding correlation function (right). The correlation strength hcorr 
is obtained from a linear ht according to Eq. [1] Since most of the statistics is at low multiplicities, the ht is mainly 
determined by the hrst points. 


their widths drj and d(p as the ALICE acceptance is approximately uniform in the selected ranges of rj 
and (p. 

2.3 Experimental procedures of the FB correlation coefficient measurement 

The present paper focuses on the study of EB correlation phenomena related to soft particle production. 
Therefore we restrict pj in 0.3 < px < E5 GeV/c, except for the study of the pj dependence presented 
in Section[5l where the pj range is 0.3 < px < 6 GeV/c. 

The correlation coefficients, ficorr. for each window pair can be calculated using two methods Ml. In 
the first method values of (ugUp), (ug), (up) and (rip) are accumulated event-by-event and then ficon- 
is determined using Eq. |2l In the second method, bcon is calculated using linear regression. The 2- 
dimensional distributions (ug. np) are obtained integrating over all selected events, then the average 
backward multiplicity is calculated for each fixed value of the forward multiplicity, and ^corr is obtained 
from a linear fit to the correlation function (see illustration in Eig. ID). Deviations from linear behavior 
may provide additional information, however, a detailed study of non-linearity in the correlation function 
is beyond the scope of this paper. 

It has been shown that the results obtained with the two methods agree within statistical uncertainty. In 
this work, results using the first method are presented. 


2.4 Corrections and systematic uncertainties 


Acceptance and tracking efficiency corrections are extracted from Monte Carlo simulations using PYTHIA6 
136] (Perugia 0 tune) and PHOJET |37, ^ as particle generators followed by a full detector response 
simulation based on GEANT3 1^ . Corrections are done to primary charged particle correlations and 
multiplicities. Correction factors obtained with these two generators are found to agree within 1% and 
the difference is neglected. Three independent correction procedures are investigated. 


In the first procedure, the correction factors for bcon- are obtained as the ratio of bcon- obtained at generator 
level (true value) to bcon after detector response simulation (measured value). In the second procedure 
the correction factors are obtained for (ugUp), (ug), (up) and (up) separately and bcon is obtained from 
the corrected moments. The third procedure takes into account approximately linear dependence of ^corr 
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Error source 

0.9 TeV 

2.76 TeV 

7 TeV 

Number of TPC space-points 

0.5-3.0 

0-0.1 

0.2-0.7 

Number of ITS space-points 

0.6-1.9 

- 

0.2-1.4 

DCA 

3.0^.0 

1.0-1.8 

O.l-l.O 

Vertex position along the beam line 

0.2-1.1 

0-1.0 

0-0.7 

^coiT correction procedure 

2.5^.0 

2.2-1.2 

1.6-2.8 

Event pile-up 

<1 

<1 

<1 

Total (%) 

3.4-1.5 

2.8^.2 

2.0-3.0 


Table 1: Sources of systematic errors of bcon measurements in Tj-windows of width Srj = 0.2, and their contribu¬ 
tions (in %). The minimal and maximal estimated values are indicated for each given source. 

on (up) when (up) varies with cuts, and each corrected value of bcon is found by extrapolation to the 
corrected value of (up). 

It was found that results of all three procedures agree within 1.6-4.2% (see Table [D, thus proving the 
robustness of bcon determination. The second procedure was chosen as the most direct and commonly 
used to produce the final corrected value of f^corr- Correction factors increase the values of bcon, obtained 
for standard cuts, by 6-10 % for analysis in T]-windows and 9-18 % for analysis in T]-0 windows and in pj 
intervals. By varying the selection cuts (vertex- , DCA- and track selection cuts), correction procedures, 
and by comparison of the high and low pile-up runs, the systematic uncertainties on bcon have been 
estimated. Adding all contributions in quadrature, the total systematic uncertainties are below 4.5% 
(4.2%, 3%) at = 0.9 (2.76, 7) TeV for the bcon analysis in T]-separated windows, and 6% for analysis 
in Tj-^ separated windows at y/? = 0.9 and 7 TeV. For the bcon analysis in pj intervals for 7 TeV, the 
systematic uncertainties are less than 8%. Statistical errors are small and within the symbol sizes for data 
points in the figures. A summary of the contributions of systematic uncertainties for bcon in t]-separated 
windows with the width Srj = 0.2 is presented in Tabled] 

3 Multiplicity correlations in windows separated in pseudorapidity 

3.1 Dependence on the gap between windows 

Fig.|5]shows the FB multiplicity correlation coefficient bcon as a function of rj^^p and for different widths 
of the rj windows (drj) in pp collisions at the three collision energies. For each y/?, bcon is found to 
decrease slowly with increasing rjgdp’ while maintaining a substantial pedestal value throughout the full 
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Fig. 5 : Forward-backward correlation strength bcon as function of rj^^p and for different windows widths 5 tj= 0.2, 
0.4, 0.6 and 0.8 in pp collisions at y/s = 0.9, 2.76 and 7 TeV. 
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bp 

Fig. 6: Correlation strength hcon- as a function of drj for = 0 in pp collisions for i/s = 0.9, 2.76 and 7 TeV. 
The MC results from PYTHIA6 Perugia 0 (solid line), Perugia 2011 (dotted line) and PHOJET (dashed line), 
calculated at generator level, are shown for comparison. The bottom panels show the ratio of bcorv between data 
and MC. The red dashed curves correspond to the model of independent particle emission from a fluctuating source 
(see text). 


T7gap range. 

3.2 Dependence on the width of windows 


The 5T]-dependence for adjacent (tjgap = 0), symmetrical windows with respect to T] = 0 is shown in 
Fig. [6l For all collision energies, the correlation coefficient increases non-linearly with S?]. This trend 
is quite well described by PYTHIA6 and PHOJET, although the agreement worsens with increasing 
This 5T]-dependence can be understood, along with other approaches J?, 25, 40], in a simple model with 
event-by-event multiplicity fluctuations and random distribution of produced particles in pseudorapidity. 
In this model, the multiplicity in an T] interval containing the fraction p of the mean multiplicity (N) in 
the full T] -acceptance is binomially distributed and its mean square is given by 


(4} = (nl}=p(l-p)(N}+p^(N^} , 

where N is the charged particle multiplicity measured in the pseudorapidity interval Y and 

(up) (ub) 5ri 

^ {N) {N) Y ■ 


(3) 

(4) 


One can connect the multiplicity fluctuations in the full rj -acceptance considered in this analysis (Y = 
1.6) with the correlation strength ^con- (see Appendix IaI): 


where 


j mod _ 
^corr 


a = 


aSrj/Y 

1 + a5r\/Y ’ 

(5) 

(A) • 

(6) 


Note that using Eq.[3]and Eq.|4]one can write the Eq.[5]also in the following form: 


j mod 1 (^f) 

^corr ^ _2 


(V) 
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(TeV) 

window width Sr] 

(np) 

^corr (^gap~i^) 

bcon (max. T]gap) 

0.9 

0.54 

1.17 

0.39 ±0.01 

0.35 ±0.01 

2.76 

0.4 

1.17 

0.44 ±0.02 

0.38 ±0.02 

7 

0.33 

1.17 

0.48 ±0.01 

0.43 ±0.01 


Table 2: Correlation strength b^on in pp collisions at ^/s = 0.9, 2.76 and 7 TeV in windows with equal mean 
multiplicity (np) and the corresponding values of Srj. Values are shown for adjacent windows (tigap=0) and for 
windows with maximal Tjg^p within Ip | < 0.8. The uncertainty on (np) is about 0.001. 

From the measured ratio of the multiplicity variance = (N^) — (A)^ in F = 1.6 to the mean value 
(JV) we obtain the value of a at ^ = 0.9, 2.76 and 7 TeV to be 2.03, 3.25 and 4.42, respectively, with 
a systematic uncertainty of about 5%. The (5Tj)-dependences calculated by Eq. [5]are shown in 

Fig.^as red dashed lines. At rj^^p = 0 the ^con-C^p) dependence is well described by this simple model. 
However, this model is not able to describe the dependence of bcon on rj^^p in Fig.[5]because it does not 
take into account the SRC contribution mentioned above. 


3.3 Dependence on the collision energy 

Figure[5]shows that the pedestal value of ^corr increases with ^/s, while the slope of the ^corrCPgap) depen¬ 
dence stays approximately constant. This indicates that the contribution of the short-range correlations 
has a very weak y^-dependence, while the long-range multiplicity correlations play a dominant role in 
pp collisions and their strength increases significantly with Note that this increase cannot be ex¬ 
plained by the increase of the mean multiplicity alone. If, at different energies, we choose window sizes 
such that the mean multiplicity stays constant the increase is still observed (see TableO. 

In the framework of the simple model described by Eq. [5] and [6] the increase of the correlation coefficient 
corresponds to the increase of the event-by-event multiplicity fluctuations with characterized by the 
ratio ay/(A). 


A strong energy dependence and rather large ficon- values were previously reported by the UA5 collabo¬ 
ration yD and recently by the ATEAS Collaboration il6ll . However, as we see in Fig.[6l the correlation 
coefficient depends in a non-linear way on the width of the pseudorapidity window. One has to take this 
fact into account when comparing the correlation strengths obtained under different experimental con¬ 
ditions. In particular, it explains the small values of bcon observed by the STAR collaboration at RHIC 
(pp, = 200 GeV) where narrow FB windows (Sr/ = 0.2) were considered, while in previous pp 
and pp experiments wider windows of a few units of pseudorapidity were used. 


4 Multiplicity correlations in windows separated in pseudorapidity and azimuth 


Multiplicity correlations are also studied in different configurations of forward and backward azimuthal 
sectors. These sectors are chosen in separated forward and backward pseudorapidity windows of width 
St] = 0.2 and d(p = n/A as shown in Fig.[3j resulting in 5 pairs with different (p-separation. 


Figs. |7] and [8] show the azimuthal dependence of bcon as a function of different T]j,gp, for 0.9 and 7 TeV, 


respectively. Data are compared to PYTHIA6 (tunes Perugia 0 and Perugia 2011), PHOJET and a para¬ 
metric string model 14111 . 


The string model fitted to our data helps to understand in a simple way the origins of the bcon behaviour. 
There are two contributions to bcon in this model. The short-range (SR) contribution originating from 
the correlation between particles produced from the decay of a single string and the long-range (ER) 
contribution arising from event-by-event fluctuations of the number of strings. The energy dependence 
of the fitted parameters demonstrates that SR parameters stay constant with ^/s while the normalized 
variance of the number of strings, the only ER parameter of the model, increases by a factor of three. 
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Fig. 7: Correlation strength bcon for separated rj-(p window pairs at = 0.9 TeV as a function of Tj separation, 
with fixed window width 5 tj= 0.2 and 5(p=n/4. The panels are for different separation distances between the two 
azimuthal sectors: (fsep = 0, 7t/4, njl, (3/4)tt and n. MC results from PYTHIA6 Perugia 0 (blue lines), Perugia 
2011 (orange dashed lines) and PHOJET (pink dashed lines) and string model 114 ill (thin green lines) are also 
shown. The bottom panels show the ratio bcon between data and MC results. 
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Fig. 8 : Correlation strength bcon for separated rj-ip window pairs at = 7 TeV as a function of 77 separation. 
The legend is the same as for Fig.|2] 


The 2-dimensional distribution of bcon as a function of and is shown in Fig. |9] for y/s = 0.9 
and 7 TeV. The qualitative behaviour of bcon resembles the results obtained for two-particle angular 
correlations: near-side peak and recoil away-side structure. The connection between the FB correlation 
and two-particle correlation function is discussed in detail in iFR 1^-431. 


The shapes of the correlation functions clearly indicate two contributions to the forward-backward mul¬ 
tiplicity correlation coefficient. The SR contribution is concentrated within a rather limited region in the 
rj-cp plane within one unit of pseudorapidity and nil in azimuth, while the FR contribution manifests 
itself as a common pedestal in the whole region of observation. 


The strength of multiplicity correlations measured in rj and rj-cp windows is compared to the results 
obtained with PYTHIA6 (tunes Perugia 0 and Perugia 2011) and PHOJET 137, 38] Monte Carlo 
generators (MC). The detailed overview of key features of these generators can be found in 1441] . Recent 
Perugia tunes for PYTHIA6 are described in i45n . 


In Fig.[T0]the comparison of bcon as a function of T]gap for 5 t] = 0.2 at ^/s = 0.9, 2.76 and 7 TeV with the 
results obtained with different MC generators is shown. All models describe the data at = 0.9 TeV 
reasonably well, while larger discrepancies are observed at 2.76 and 7 TeV, with PYTHIA giving a 
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Fig. 9: 2D representation of boon- at (a) ^/s — 0.9 TeV and (b) at ^/s = 1 TeV for separated r\-<p window pairs with 
5 tj= 0.2 and d(p=7t/4. To improve visibility, the point (tisep’?’,sep)=(0’0) and thus i7corr=l is limited to the level of 
the maximum value in adjacent bins. 
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Fig. 10: Correlation sttength bcorr as a function of rj^ap in pp collisions for data taken from Fig.|5]and compared to 
MC generators PYTHIA Perugia 0 (solid line), Perugia 2011 (dotted line) and PHOJET (dashed line) for = 0.9, 
2.76 and 7 TeV collision energies, windows width Srj is 0.2. The bottom panels show the ratio of the data to MC. 

better description of the data than PHOJET. Qualitatively similar conclusions can be drawn from the 
comparison of the 5T]-dependence in experimental data and MC as shown in Eig.0 

Note that PYTHIA also describes the correlations in rj-q) windows reasonably well, see Eigures |7] and 
[H while PHOJET gives a good description only for ^/s = 0.9 TeV and significantly underestimates the 
data at 7 TeV. 

The difference between the experimental data and the results obtained with MC generators is more visible 
in Eig.[TT] which compares the measured ratio of bcon at = 2.76 and 7 TeV with respect to 0.9 TeV 
as a function of rj^^p to MC calculations. The measured ratios show an increasing trend as a function of 
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Fig. 11: Ratio of /jcorr at 2.76 (blue squares) and 7 TeV (red circles) with respect to 0.9 TeV vs. Pgap- The 
calculations from MC generators are also shown: PYTHIA Perugia 0 (solid line), Perugia 2011 (dotted line) and 
PHOJET (dashed line), for 2.76 TeV (blue lines) and 7 TeV (red lines). 


T]gap, while PYTHIA and PHOJET underestimate the ratios and exhibit a flatter T]gap dependence. 

It is important to note that, in the framework of PYTHIA, the observed ER part of bma: (the pedestal 
in Eig. |2l is dominated by multiple parton-parton interactions (MPI). This supports earlier results 
in which the EB correlations in pp collisions were studied by MC simulations with recent tunes of the 
PYTHIA6 at ^/s = 0.9 TeV. Hence, the observed dependence of bcon on collision energy and on different 
configurations of rapidity and azimuthal windows adds new constraints on phenomenological models for 
multi-particle production. 


5 Dependence of FB multiplicity correlation strength on the choice of pj intervals 


The behaviour of EB multiplicity correlation strength was also studied as a function of pT of registered 
particles. These studies were motivated by a recent paper by the ATEAS collaboration ifltill . which 
reported a decrease in the multiplicity correlation strength with increasing p™". However, as we have 
observed in Section 13.21 there is a strong non-linear dependence of b^ou on the size of pseudorapidity 
windows and, hence, on the mean multiplicity (uch) in the window (see Eq. |4l|5l and Eig.lQ. In order 
to demonstrate that the strong p™" dependence is not a trivial multiplicity dependence, in our analysis 
we use pt intervals with the same (uch). To this end, the correlation strength bmu is studied for five px 
intervals within 0.3 < px < 6 GeV/c at y/s = 7 TeV: 0.3-0.4, 0.4-0.52, 0.52-0.7, 0.7-1.03 and 1.03-6.0 
(GeV/c). In each px interval, the corrected mean multiplicity (up) = 0.157 with a systematic uncertainty 
about 2%. Correlations are studied in rj and rj-q) EB-windows configurations. Note that in case of 
windows chosen symmetrically with respect to T] = 0 the definition of ^con- given by Q coincides with 
the correlation coefficient Ppg used in the ATEAS analysis. 


Eig.[l2] shows bcon as a function of pj“ for rj^^p = 0 and 1.2. Systematic uncertainties are shown as 
rectangles, statistical uncertainties are negligible. We find that bcon increases with pj“ for both values 
of figap, in contrast to the results reported in iltill . This result can be understood if one takes into ac¬ 
count that the multiplicity fluctuations in a given window are closely connected with the two-particle 
correlation strength 10, S]. In the simple model with the event-by-event multiplicity fluctuations and 
random distribution of produced particles in pseudorapidity, discussed in Section [3^ Eq.|7]allows us to 
discuss the observed dependence of the correlation coefficient bcon on the px-binnings for the case of 
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Fig. 12: Correlation strength boon- at = 7 TeV for separated pseudorapidity window pairs, measured in pj 
intervals with same (wch), as a function of p™" for each interval. Values are shown for Pgap=0, 1.2 with 5 tj=0.2. 
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Fig. 14: Correlation strength bcon- at ^/s = 1 TeV for separated r\-<p windows in different pi intervals with same 
(wch)- Five values are shown as a function of rj^gp. Windows of dr] =0.2 and 5<p=n/4. Top and bottom panels 
and contain the same experimental data, lines correspond to PYTHIA6 Perugia 2011 (top) and to PHOJET (bot¬ 
tom). 


Vgap = 0 (Fig.fT^. One sees that the imposed condition {np)=const eliminates the dependence of bcon 
on the multiplicity. The ratio 1 /decreases and increases with increasing p™ 

As mentioned above, in the approach used in J3] the dependence of the correlation strength on the 
of charged particles was studied without cuts on which leads to a decrease of the correlation 
^corr with increasing p™". This result can also be illustrated with the help of Eq. |7J In this case (up) 
decreases with increasing pj “ and (tip) /increases (approaching the Poisson limit = (rip)) leading 
to the decrease of bf°^. Thus, the difference of the results in these two approaches can be qualitatively 
understood using Eq. |7J 

Eig. [I3] shows Z^corr as function of for different pj intervals. Eig.[T3](a) compares data to PYTHIA6 
tune Perugia 2011. The general trend of bcon increasing with higher pj“ for all is reproduced by 
this tune, with small quantitative deviations. Pig.[l3](b) shows the same data in comparison to PHOJET. 
This generator does not describe the data well: PHOJET results are almost independent of pjand only 
grow significantly for the px range 1.03-6.00 (GeV/c). Since experimental data was used to determine 
the Px intervals with the same mean multiplicity, the values of mean multiplicities may vary slightly in 
case of the MC samples for the same px intervals. Deviations from the mean value are within 4% for 
PYTHIA6 Perugia 0 and 12% for PHOJET. 


The analysis of 6 coit is also performed in rj-rp separated windows in different px intervals with the 
same mean multiplicity (for pp collisions at ^ = 1 TeV) in 8x8 rj-rp windows. Results are shown in 
Eig.[l4]and compared to PYTHIA6 and PHOJET calculations. In addition to the conclusions that were 
drawn above from the correlations between T]-separated windows, some new details are revealed. In 
particular, one observes that the PHOJET discrepancy with the data is especially dramatic at (p^^p=nl2, 
where PHOJET shows no dependence of bcon on the px range. It was shown already in i47n that PHOJET 
has difficulties in description of underlying event measurements. 

Eig.[l4] shows that for higher px intervals a near-side peak appears (see panels for <Psep = 0 and 7r/4), 
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at the same time the bcoii in the flat region at ri^ep > 1 increases with px for all <Psep values (compare 
panels for <Psep = 7t/2, 37r/4 and n). It should be emphasized that the value of the pedestal (the common 
constant component in all panels) increases with px- 

In near- and away-side azimuthal regions the increase of bcoii with pjp“ can be explained by an enhanced 
number of back-to-back decays and jets. The general rise of bma: can be related to the increase of the 
variance in Eq.[6l discussed in the framework of the simple model in Section 

6 Conclusion 

The strengths of forward-backward (EB) multiplicity correlations have been measured in minimum bias 
pp collisions at ^/s = 0.9, 2.76 and 7 TeV using multiplicities determined in two separated pseudorapidity 
windows separated by a variable gap, T]g^p, of up to 1.2 units. The dependences of the correlation 
coefficient bcoa: on the collision energy, the width and the position of pseudorapidity windows have 
been investigated. Eor the first time, the analysis has been also applied for various configurations of the 
azimuthal sectors selected within these pseudorapidity windows in events at y/s = 0.9 and 7 TeV. 

A considerable increase of the EB correlation strength with the growth of the collision energy from 
^/s = 0.9 to 7 TeV is observed. It is shown that this cannot be explained by the increase of the mean 
multiplicity alone. The correlation strength grows with the width of pseudorapidity windows, while it 
decreases slightly with increasing pseudorapidity gap between the windows. It is shown that there is a 
strong non-linear dependence of the correlation strength on the width of the pseudorapidity windows and 
hence on the mean multiplicity value. 

Measurements of the correlation strength for various configurations of azimuthal sectors enable the dis¬ 
tinction of two contributions: short-range (SR) and long-range (ER) correlations. A weak dependence on 
the collision energy is observed for the SR component while the ER component has a strong dependence. 
Eor T]-gaps larger than one unit of pseudorapidity and 71/2 in azimuth the ER contribution dominates. This 
contribution forms a pedestal value (the common constant component) of b^ou increasing with collision 
energy. 

Moreover, pseudorapidity and pseudorapidity-azimuthal distributions of bcon have been obtained in pp 
events at ^/s = 7 TeV for various particle transverse momentum intervals. It is found that the EB corre¬ 
lation strength increases with the transverse momentum if px-intervals with the same mean multiplicity 
are chosen. 

The measurements have been compared to calculations using the PYTHIA and PHOJET MC event gen¬ 
erators. These generators are able to describe the general trends of b^ou as a function of drj, rj^^p and (p^gp 
and its dependence on the collision energy. In px-dependent analysis of b^on, PYTHIA describes data 
reasonably well, while PHOJET fails to describe bcon in azimuthal sectors. The observed dependences 
of bcon add new constraints on phenomenological models. In particular the transition between soft and 
hard processes in pp collisions can be investigated in detail using the pj dependence of azimuthal and 
pseudorapidity distributions of forward-backward multiplicity correlation strength /?corr- 
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A A model with random uniform distribution of produced particles in pseudorapidity 


In a simple model with event-by-event multiplicity fluctuations and random uniform distribution of pro¬ 
duced particles in pseudorapidity the probability to observe n^ particles in some subinterval St] from 
the total number of N charged particles produced in the whole pseudorapidity interval Y is given by the 
binomial distribution: 

Pv(%)=C>"^(l-p)'^-"F , (A.l) 

with (np}jy = pN and = p{\ — p)N -|- p^N'^, where p = 5r]/Y. (We consider the case of symmetric 
windows Srjp = Srjg = St].) Averaging then over events with different values ofV, 

P(np) = £p(V)Pv(«f) , (A.2) 

N 


we have 


and hence 


{np) = Y^P{np)np = ££P(V)PvK)% = Y,PiN)pN = p (V) 

Mp ”f ''V N 

(N) (N) Y • 


In the same way we find 

(4) = («b) =P(1 -P){N) + , 

{{np + npf) = 2p{\- 2p) {N) + {2pf{N^) . 


One can rewrite (IA.4I) - (IA.6I) also as 

~ + np) _ (^1 - (np) _ (J^-{N} _ 

(up + nj^y {np}^ (V)2 “ ’ 


(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 
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since the so-called robust variance /?n is the same for any subinterval of Y in the case of the independent 
homogeneous distribution of the particles along Y i43n . 

Using the presentation for the covariance 


(Uf^b) “ (%)W) = ~ ’ 


we can write for the correlation coefficient in a model-independent way: 


^cniT — 


-^np+ng '^tip 




2o}_ 


Then combining (IA.7I) and (IA.9I) we find 


j mod _ 
^corr 


(»F)f?N 

1 + 


Using (IA.4I) we can write (lA.lOl) also as 


, mod ^ (^Sri/Y 

l + a5T7/T 


where 


a = (A^)/?n 



Substituting the expression 


/?N = 



from (IA.7I) into (lA.lOb one finds another presentation for : 


j mod 1 (^f) 

^con' ^ _2 


(A.8) 

(A.9) 

(A. 10) 

(A. 11) 

(A. 12) 

(A. 13) 

(A. 14) 
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